The extended Lawrence-Doniach model: the temperature evolution of the in-plane 

magnetic field anisotropy 



(N 

o 

(N 

> 
O 



o 
o 

u . 
=5 

a- 



-a 

o 
o 



> 

(N 



(N 



X 



M. D. Croitoru 1 and A. I. Buzdin 1 - 2 
1 Umversite Bordeaux I, LOMA, UMR 5798, F-33400 Talence, France and 
2 Institut Unwersitaire de France, Paris 
(Dated: November 28, 2012) 

Using the quasi-classical formalism we provide the description of the temperature and field- 
direction dependence of the in-plane upper critical field in layered superconductors, taking into 
account the interlayer Josephson coupling and the paramagnetic spin splitting. We generalize the 
Lawrence-Doniach model for the case of high magnetic fields and show that the re-entrant supercon- 
ductivity is naturally described by our formalism when neglecting the Pauli pair breaking effect. We 
demonstrate that in layered superconductors the in-plane anisotropy of the onset of superconduc- 
tivity exhibits four different temperature regimes: from the Ginzburg-Landau type in the vicinity of 
the critical temperature T c o with anisotropies of coherence lengths, up to the FFLO type induced 
by the strong interference between the modulation vector and the orbital effect. Our results are in 
agreement with the experimental measurements of the field-angle dependence of the superconducting 
onset temperature of the organic compound (TMTSF^CICU. 

PACS numbers: 



I. INTRODUCTION 

Since the discovery of superconductivity in the first 
layered compound^ there have been found many types 
of superconductors consisting of alternating conducting 
and insulating layers. Examples include the high-T c 
cuprates^ layered ruthenates^ the iron pnictides and 
oxypnictides^— , graphite intercalation compounds^— 
crystalline organic metals* 1 ^ - — the various types of ar- 



tificial multi-layers 



14-16 



etc. Amongst them, layered 



organic metals arc distinctive for a number of reasons. 
Most of them exhibit profound reduced dimensionality 
reflected in the very strong charge-transfer anisotropy. 
The interplay between electronic correlations and en- 
hanced dimensionality effects leads to a broad range of 
physical properties observed in these materials. More- 
over, organic metals are often available in highly clean 
single crystals that enables one to perform detailed band- 
structure measurements and to study mechanisms of su- 
perconductivity in quasi-low-dimensional electronic sys- 
tems. Finally one of the most prominent property of or- 
ganic layered superconductors is their robustness against 
high magnetic fields applied parallel to the conduction 
layers. Commonly known examples include Bechgaard 
salt superconductors (TMTSF^X, where anion X is PF@, 
CIO4 etc. A very large upper critical fields, which exceed 
the Pauli paramagnetic limit, for a magnetic field aligned 
parallel to their conducting layers were reported* 1 ^— In 
the compound (TMTSF) 2 PF 6 H c2 = 90 kOe^ which 
is more than 4 times larger than Hp ~ 22 kOe and an 
enhancement of almost two times over Hp ~ 27 kOe 
is observed in the compound (TMTSF) 2 C10 4 , H c2 ^ 
50 kOe» 

In magnetic field the superconductivity in usual type 
II superconductors is suppressed due to the diamagnetic 
currents and the Pauli pair breaking effect for singlet 
pairing. In layered conductors the spatial orbital motion 



of electrons is mostly restricted to the conducting planes, 
when charge carrier hopping between adjacent layers is 
small, and the magnetic field applied precisely parallel to 
the conducting planes weakly affects the orbital motion of 
electrons. Hence the orbital depairing is largely avoided 
(there is no magnetic flux inside the 2D Cooper pairs lo- 
cated in planes in such situation). Moreover, when the 
interlayer coherence length in a quasi-lD superconductor 
is comparable to the interlayer distance the field-induced 
quasi-2D (3D) — » 2D dimensional crossover occurs in a 
high magnetic field, restoring the bare critical tempera- 
ture, T c0 mi 

Various theories based on different pairing symme- 
tries predicting the existence of high-field superconduct- 
ing state have been proposed previously. Among them, 
a phase transition to an inhomogeneous FFLO (Fuldc- 
Ferell-Larkin-Ovchinnikov) phase for T < T* ~ 0.56T c0 
or H > H* ~ 1.06T c0 /hb, in which the the sin- 
glet superconducting ground state is characterized by 
the spatially modulated order parameter and the spin- 
polarization. Therefore superconducting state can be 
stable beyond the field set by the Pauli paramagnetic 
limit, [inHp = Aq/\/2, where Ao is the superconduct- 
ing gap at T = . 23 i 24 Conditions for the stabilization of 
the FFLO phase are rather stringent^ namely (i) the 
orbital pair breaking effect should be sufficiently weaker 
than the Pauli paramagnetic limit, the Maki parameter 
Q'm = V2H C 2/ Hp > 1.8; (ii) the system should be in 
a clean limiti2£r— A growing body of experimental evi- 
dence for the FFLO phase reported from a various mea- 
surement techniques supports this scenario^r— An al- 
ternative to the FFLO phase is a triplet pairing state, 
when the Pauli spin-splitting destructive mechanism is 
absent. Within this pairing symmetry, as was shown 
by Lebed , 22 ' 46 the superconducting state is always sta- 
ble at low temperatures and exhibits a strong re-entrant 
behavior in high magnetic field. So far the re-entrant su- 



perconducting phase has not been experimentally iden- 
tified, at least it is difficult to make more than a tenta- 
tive judgement.— Nevertheless it can reveal itself in a 
number of nontrivial effects in singlet-paired organic ma- 
terials in high magnetic fields i 52 ' 53 It was shown that it 
can appear in a hidden form and be responsible for an 
increase of the superconducting transition temperature 
in a magnetic field if the orbital effects of an electron 
motion are stronger than the Pauli spin-splitting effects 
(Paramagnetic intrinsic Meissner effect).— 

Hitherto there is no experiment which unequivocally 
answer on the lingering question concerning the su- 
perconducting pairing symmetry in (TMTSF^X com- 
pounds. Previously it was reported that the Knight 
shift in (TMTSF) 2 PF 6 conductor does not change 
at transition temperature supporting the triplet sce- 
nario of pairing.— However, later experiment with 
(TMTSF)2C104 conductor at low-field regime have re- 
vealed a clear change of the Knight shift at the super- 
conducting transition making possible consideration of 
the singlet scenario of pairing in such structures^ In 
the high-field regime the Knight shift is quite weak. On 
the other hand, as shown in Ref.— a small fraction of 
the triplet pairing in the singlet paired superconductor 
strongly enhances the upper critical field and the triplet 
component of the order parameter is always generated in 
singlet superconductors due to the Pauli paramagnetic 
spin-splitting effects . 48 ' 49 

In this paper we extend results presented in our previ- 
ous Letter,— and investigate the in-plane magnetic field- 
angle dependence of the onset of superconductivity in 
layered conductors in the conventional and the FFLO 
modulated phases. For this purpose we provide the qua- 
siclassical description of the anisotropy of the in-plane 
critical field in layered superconductors and generalize 
the Lowcrcncc-Doniach model for the case of high mag- 
netic fields. 

The layout of our paper is as follows. In Scc|Hl wc 
outline our model based on the quasi-classical formalism 
for layered superconducting samples. In Sec. IIIII we de- 
rive the generalized Lowerence-Doniach equation. In sec. 
IIVI we extend this model to the extremely high magnetic 
fields. In Sec. |V] we focus on the in-plane anisotropy of 
the upper critical field for layered superconductors when 
only orbital motion is included in the model, and then we 
investigate the in-plane anisotropy of H C 2 when both or- 
bital and paramagnetic depairing are accounted for. Fi- 
nally, a short summary is given, where we emphasize the 
significance of the obtained results for the interpretation 
of experiments with layered superconductors. 



II. GENERAL SETTINGS 

We consider a system consisting of layers with good 
conductivity in xy-planc stacked along the z-axis [see 




FIG. 1: Scheme of the quasi-2D layered metal. 



Fig- ID- The single-electron spectrum is taken as follows 
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where e(p z ) = It cos (p z d) with d - the interlayer dis- 
tance. We assume that the coupling between layers is 
small [see Fig. [5], i.e. t <C T c0 , but sufficiently large to 
make the mean field treatment justified, T^ /Ep <§; 
Here T c q is the critical temperature of the system at 
H = 0. In purely 2D samples, phase fluctuations de- 
stroy the long-range order, however as shown in Ref^, 
even a very small value of hopping leads to restoration of 
superconducting order. 

We choose the magnetic field to be parallel to the con- 
ducting planes and with a gauge for which the vector 
potential A = H x r [r = (x, y, 0) is a coordinate in .De- 
plane], i.e. A z = —xHsina + yH cos a, where a is the 
angle between the applied field, with amplitude H, and 
x-axis. Assuming that the vector potential varies slowly 
at the interlayer distances (this assumption means that 
we neglect the diamagnetic screening currents and take 
the magnetic field as uniform and given by the external 
field, H ) , and taking into account that the system is near 
the second-order phase transition, we can employ the lin- 
earized Eilenberger equation for a layered superconduc- 
tor in the presence of the parallel magnetic field (in the 
momentum representation with respect to the coordinate 
z)2L: 



n„+fi) U(n,r,p z ,k z ) = [A{r,k z ) 
(f u (n,r, p z ,k z )) 



2t 



sign(w„). (2) 



Here 



~ h k 

n= -v F .V + 2itsm(ps«2)sin(Q.r — ^-d), (3) 

where Q = (-KdH/4> )[— sina, (mjmj) 1 ' 2 cos a, 0] with 
4>o = irhc/e, h = is the Zeeman energy, vp = u^n 

is the in-plane Fermi velocity, r is the impurity scatter- 
ing time, and il„ = ui n — ih + sign(a;„)/2r. The order 
parameter is defined self-consistently as 



(r, k z ) = 2nm ^ (f u (n, r,p z , k z )) 



(4) 
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FIG. 2: The Fermi surface of the layered metal in the form of 
a corrugated cylinder. 



To obtain the convergent iterative scheme we need to re- 
quire that A(r) 3> frv p. A(r) /2ttT c o, that implies that 
characteristic scale of the order parameter variations 
should be much larger than £q = hvp /2ttT c q. After the 
completion of the k-th iteration we obtain 

(n, r, Pz ,k z ) = J2 ^TT fii A(r, k z ). (9) 



Taking into account the averaging procedure over mo- 
mentum p z , and hence omitting the terms with even 
powers of sin (p z d) , then retaining terms up to the sec- 
ond order in fivi7.VA(r)/27rT c o, and making use of the 
self-consistency relation Eq. (j4|), we obtain the ex- 
tended Lowerence-Doniach equation (MLD equation) in 
the isotropic case 



where A is the pairing constant and the brackets denote 
averaging over p z and n, 



2tt 

d dp z f da 



2tt J 2tt 
o 



(...). 



(5) 



We assume that the temperature unit is so chosen that 
the Boltzmann constant kp = 1. 

Here we considered a layered superconductor in the 
clean limit, meaning that the in-plane mean free path 
is much larger than the corresponding intra-plane coher- 
ence length, £q = hvp / (2ttT c q). Therefore the linearized 
Eilenbcrgcr equation for the anomalous Green function 
/ w (n, r,p z ,k z ) describing layered superconducting sys- 
tems acquires the form 



0„ + n f u (n,r, Pz ,k z ) = A(r,k z ) (6) 



with Q n = uj n — ih sign(w ra ) from now on. 



III. A LAYERED SUPERCONDUCTOR IN A 
PARALLEL MAGNETIC FIELD 

The upper critical field corresponds to the highest 
value of H, for which the solution of Eqs. (j4|) and ((6]) 
exists. To start with, we consider Eq. (j6]) and write it in 
the form 

f u (n,r,p z ,k z ) = A ^ kz ^ - J-ILf u (n,T,p z ,k z ) (7) 

convenient for the subsequent derivation of iterative pro- 
cedure. Using this equation we construct the following 
iterative scheme 



(n,r,p„k,) 



A(r,fc 2 



i 



n/( fc )(n,r,p 2 ,fc z 



(8) 



A (r, k z ) In 



A (r, k z 



1 



n 



h=£0 
MLD 



A(r,k z ) (10) 



where T c q is the critical temperature in the absence of 
coupling between adjacent layers, t, and of the magnetic 
field, and 



n 



MLD 



n n 

8 



^3 [1 - cos (2Q.r-M)] 



[1 - 7cos(2Q.r-fc z d)] (11) 



The anisotropic case one can obtain simply 
by the following substitutions h 2 v F W 2 A (r) — > 
2e(V)A(r) = 2h 2 {(v 2 Fx ) d 2 + (v 2 Fy ) d 2 } A (r) and 

H 2 (v F Q) 2 -> 2e(Q), where 

e(Q)^h 2 {{v 2 Fx )Q 2 x + {v 2 Fy )Q 2 y }. (12) 

Introducing the temperature T c p, as the superconducting 
onset temperature in the pure Pauli limit determined by 
the expression 



hi 



I c0 



i 



i 

{ 1,, 



(13) 



and the use of the identities 2nTJ2n=o^n = 
-$( 2 ) (/i) /8^ 2 T 2 and 27rT^ =o r!- 5 
-$( 4 ) (h) /3847T 2 T 2 gives rise to (for details see 
Appendix [AJ 



A (r, k z )P = - 



$( 2 ) (h) 

8ir 2 T 2 



< - t 2 [1 - cos (2Q.r-fc,d)] } A (r, k z ) 



$< 4 > (h) t 2 e (Q) 
3847r 4 T c 4 p 4 



[1 - 7cos (2Q.r-fc z d)] A (r, k z ) , 



(14) 



4 



where P = (T c - T cP ) /AT C , $( fe) (h) = 
[V> (fc) (l/2 + i/i) + V' (fc) {1/2- ih)] /2 with V (fe) (^) = 
(2:) /dz k and ^ (z) is the digamma func- 
tion. If we can neglect the Zeeman effect, 
h = 0, than Q„ — > u„, and making use of the 



identities 2irT Y^= 
2vrTEr=o^ 5 ' 



7C (3) /4ir 2 T 2 and 
31C (5) /16tt 2 T 2 reduces Eq. dHJ) 



to 



A (r, fc z ) In ^ - 47r 2y2 Q 

< ^—2 - t 2 [1 - cos (2Q.r-M)] }• A (r, 



31C(5)f 2 £ (Q) 
16^ T 4 4 



[1 - 7 cos (2Q.r-fed)] A (r, fc z ) , 

(15) 
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where T c q is the superconducting critical temperature in 
the absence of coupling between adjacent layers, t, and 
in the absence of the magnetic field, described by the 
vector Q. As it is seen the MLD equation contains the 
term, proportional to (vpQ) 2 t 2 7 which is absent in the 
standard Lowcrence-Doniach equation. As it will be seen 
later this term represents unusual orbital contribution 
responsible for the re-entrant superconducting phase at 
high magnetic fields^ 2 - Let us consider several limiting 
cases. 



A. Regime H <C 



First, let us consider the case of a small magnetic field. 
When hvpQ <C T c0 , we can retain only terms up to the 
second order in {HvpQ) /T c q or/and t/T C Q. Then after 
neglecting the last term in the MLD, because it is much 
smaller than other terms, Eq. (fTOj) reduces to the stan- 
dard Lowerence-Doniach equation 



A(r,fc 2 )P = ^T c £lgU(r,fc 2 



^2 sin 2 (Q.r-^d 



In the continuous limit, d — > 0, d <C £q (-0 with £g- 
- the inter-plane coherence length, Eq. (|16[) transforms 
into the Ginzburg-Landau equation for an anisotropic su- 
perconductor. If the order parameter is homogeneous 
along the 2-axis we can set k z = 0. If Qr^Ql <C 1, 
where I 



or H -C —A — d>o, 
tcristic magnetic length with the characteristic magnetic 



y^h/mu>H is the charac- 



frequency ujh defined as ujh 
can be further simplified 



m.„ LU (ha > 



Eq. (H6 



PA (r) 



lxd 2 x + jyd 2 - "f z 



2Qr 



A(r) = 0, (17) 



FIG. 3: Scheme of the H — T phase diagram for layered su- 
perconductors, when external magnetic field, H, is applied 
parallel to the layers, t <C T c o, and the paramagnetic effects 
are vanished. 



■ 2 ,)/32n 2 T 2 



where a 

-h 2 & 2 ) {h) (v F 

d 2 t 2 ^ (h) /16n 2 T 2 P . If h = 
7x,„ = Ph 2 (v 2 Fx J /2T% = 
where (3 = 7C(3)/8tt 2 . 
is ujh = \H^ — irH, or using 
l x h z = (v 2 Fx J/2d 2 t 2 = v 2 p/4d 2 t 2 



[T c — T c q) /T c o, 

lz = 

we may write 
2 4/4T 2 , 72 = Pd 2 i 2 /T 2 a , 
The cyclotron frequency 
the relation 

_ 2dt 2tt 



IS UJH 



-H. 



After performing scaling of the variable y' = y 'm y / 'm x y , 
the anisotropic model with effective masses can be 
reduced to the isotropic one in the renormalized mag- 

1 where 



netic field H 

X2 



^sin 2 (0) + ^ cos 2 {■&), 

ti 2 /2m x ^ y = jx.yTco- Finally, the angle-resolved highest 
magnetic field, at which superconductivity can nucleate 
in a sample is given by 



Hc2 (0,T)|, 



H «(?)\ 



sin 2 {■&) + ^ cos 2 (1?) 



(18) 



Here for the negligible Zeeman effect, which breaks apart 
the paired electrons if they are in a spin-singlet state, 
h = 0, 



H, 



h=0 



h 2 d t 2tt 1 T c0 



(19) 



while for h 7^ 



H: 



h^0 



8ttT, 



m r T. 



where xi : H <C 



Ahdt V 2$( 2 ) (h) 2tt 



1 - 



2k 

T, 



cP 



(20) 



5 



-Tco 



B. The crossover regime: 



To study the anisotropy of the upper critical field, 
when its amplitude is in the range, t <C HvfQ <C T c q, 



-t < P < 



<t>a ± 
■nhdvp 



T c q, we employ the extended 



Lowerence-Doniach equation Eq. (fTTJ|) and choose the 
solution in the form 

A(r) = A + A 2 cos(2Q.r). (21) 

By substitution it in Eq. (TTU)) , wc obtain the following 
system of coupled equations 



A P = nT cP ^ 



1 

E 



g(Q) 

1 



i 2 A 
e(Q) 



7 



80S 



t 2 A 2 , 



(22) 



2.0 



?r 1.5 
-s; 



1.0 
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45 



0.0 



/7=3.2 



^2.55 



f= 1.7 



FFLO 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



and 
A 2 P 



5A 2 \ e(Q)i 2 



8^ 



(23) 
when 



In the situation |P| < |$( 2 ) (h)\ e (Q) /8tt 2 T- 
taking into account that Ao 3> A 2 , from Eq. (f2"3")l we 
can obtain A 2 = t 2 Ao/e(Q). Substituting it into Eq. 
(|22[) and retaining only terms up to the second order in 
(t/T c0 ) leads to 



P 



7rT c p2j 



4fi5 2fi3 £ (Q) 



(24) 



FIG. 4: hqvF /kBTco and HQvf /fcsTcO versus T c p/T c o for sev- 
eral values of 77. 



J. Lowerence-Doniach regime 
t4>o/nhdvF <C P <C V~tTco4>o /nfidvp 

If the amplitude of the external magnetic field satisfies 
the condition t <C SwfQ -C VtTcO, or tcjjo/TrhdvF -C P <C 
\/tT C Q(f>o /irhdvp , we can neglect the first term in Eq. (f2"u]l 
and obtain Eq. (|18[) . as the expression for the upper 
critical field with 



P, 



fj=0 
c2 



2m T 



7C(3)t 2 
8?r 3 MT c0 V T c - T f 



(29) 



$( 2 ) (fe) 

8^; 



t 2 -e(Q)t 



2 $ (4) w 



t 4 $( 2 ) (ft) 



ISSBtt 2 ^ e (Q) 167r 2 T c 2 F 

(25) 

If the Zccman effect of the applied field is absent, /i = 
0, we have to make the following substitution, fi„ — > 
w„ and T c p — > T c q. After introducing the temperature, 
T c t, accounting for the coupling between adjacent layers 



via expression ln(T c0 /T ct ) = t 2 itT ct Y^ n - 
acquires the form 



,-3 



Eq. 



Tct 4" 4^5 



where 



e(Q) 



fl 2 V 2 F TT 2 d 2 H 2 

2 <t> 2 



sin 2 (tf) + 01* cos 2 (#) 



or using the definition (v Fxy ) = 



e(Q) 



ft 2 T rn ^ 2 d 2 P 2 



/3 <t>, 



m , 



sin 2 (7?) + — cos 2 (7?) 



(26) 



(27) 



(28) 



Eq. (|26[) is the transcendental equation to determine 
H C 2 (1?) 7 1 ) for a layered system with interlaycr coupling 
t. Let us consider two limiting situations. 



where yen : tifio/nhdvF <§C P *C s/tT c o<t>o /irhdvF , when 
the Zeeman effect is negligible. 



2- ^ff^^ 7 V«k«P<^T c o 



If the field is such that y/tT c o <C HvfQ <C P c o, or 
^t^VW^ « H « iM^^cO the expression for the 
upper critical field, P c2 (i?,T), can be obtained from Eq. 



(|26[) by neglecting the second term. Then again we obtain 
Eq. (fT8|). as the expression for the upper critical field 
with 



P 



h=Q (K 

c2 \2 



31 C (5) dt 



28C(3)P c0 2m x (T c -T ct ) 



r, 2 



tu- 



rn 

This regime describes the beginning of the reentrant su- 
perconductivity regime i 22 ' 46 



IV. GENERAL CASE FOR H » 



: (j)Q 



To study the anisotropy of the upper critical field, 
when its amplitude satisfies P 3> Khdv F 0o wc need 



6 



to reconsider the solution of the Eilcnbcrger equation, 
Eq. Since the magnetic field induced poten- 

tial has the form V (r) = tsm(p z d) [e lQ - r - e'^] = 
lit sin (p z d) sin (Q.r), i.e. it is periodic in real space, the 
solution of Eq. (|6]) can be written without any loss of 
generality as^ 

L (n p ,r, Pz ) = e^J2 e " nQr fm K,n p ,p z ) , (31) 



where we took into account the possibility of the FFLO 
phase formation in this field regime. Because of the form 
for f u (n p , r,p z ) of Eq. (|3"Tj) one can write A (r) as 



A(r) = e iq r ^V 2mQ r A 2li 



(32) 



From symmetry considerations it follows that A_ 2m = 
A 2m . Substituting Eqs. ((31]) and (J32J) back into Eq. © 
one gets^S 

Ln(q)/o + i/-i-i/i = A , (33) 
L n (q ± Q) /±i±t/o T tf± 2 = 0, (34) 

Mq±2Q)/±2±t/±i=F*/±3 = A± 2 , ( 35 ) 
L„(q±3Q)/ ±3 ±t/±a = J (36) 

where /,„ = f m (w„, n,p z ), L n (s) = fi„ + ihv F s/2 
and t = tsm(p z d). Here we took into account that 
A+ 

(2m+i) = 0. When deriving this set of coupled equa- 



tions we accounted for t <C hvpQ, or 



00 



t < H. This 



limit allowed us to retain only Ao and A±2, or /o, f±i, 
f±2 harmonics, because we adopt a second-order approx- 
imation in the small parameter t/T c o to the solution of 
Eq. (0, t <C T c o). Actually, if the applied field is such 
that T c o < HvfQ, then it would be sufficient to retain 
only Ao, or /o, f±i harmonics. 

Making use of the self-consistency relation the solution 
of the system of coupled equations (|33l - l3"6")l can be given 
in the form (for details see Appendix [Bj 

A [P + t 2 a] =; 2 5>±A ±2 , (37) 



A +2 [P + t 2 b + +S+] =t 2 c+A , 



A_ 2 \P + t 2 b_ + S-] = i 2 c_A 



(38) 

-2 |/~ -t-t o- -i- o-j = t c_z-io (39) 
where the following notations are introduced: 

a = irT r„(q,q,£Q)| T=Tep) (40) 



n,£=± 



& ± =7rT^] T„(q±2Q,q±2Q,q±2Q+£Q)| 



C± =^T^T„(q,q±Q,q±2Q)| 

n 



T=T C . 



L n (q) L„(q±2Q) 



T=T cP ' 

(41) 
(42) 

(43) 



T=T C , 



with T n (q, p, k) = (L- 1 (q) L- 1 (p) i,- 1 (k)> /2. The so- 
lution of the system (|37 |) - (|38l) is found from 



P + t 2 b- 

-t 2 c_ 




5_ -t 2 







P + i 2 a 

-t 2 c+ P + t 2 b+ 



t 2 c+ 



0. (44) 



For T > T*, when q = 0, A +2 = A_ 2 , which makes it 
possible to write the solution in the form 



with 
S ± (Q) 



T c = T cP [1 - AS ± (Q)] 



(a + 6±)t 2 + 5^ 



(45) 



2 



4c±^c±. (46) 



If VtT c o <C S^fQ then it further simplifies, S ± (Q) = 
at 2 . In Eq. (|45|) those values of ± are chosen that maxi- 
mize the critical temperature. In general case, if H < ff* 
then within a second-order approximation in the small 
parameter £/T c o, A±2 reads as 



A 



■Ao 



(47) 



and the solution (|45j) system (|37fl38f) simplifies to (for 
details see Appendix [Cjl 



_ 1+ i(hv F QY 



<>7 



(48) 



In the absence of the Zccman effect 



Tot 



7T 2 T 2 



31C(5 ) (fapQ) 2 , 7C(3) 
128 



7T 2 T 2 



(hvFQ) 



(49) 

which is the same as Eq. (|26|) . Thus, within the expan- 
sion model (|31[) we obtained the upper critical field versus 
the superconducting onset temperature. This equation 
naturally describes the crossover between two regimes: 
the Lowerence-Doniach phase and the beginning of the 
Lebed re-entrant phase. 



3. Regime of high magnetic fields (jHf 3> w ^2v F < ^ >0 ) 

In the absence of the Zeeman effect, when studying 
the anisotropy of the upper critical field, such as T c q <C 
hvpQ, the second harmonics in the expansion Eq. (|31[) 
can be neglected, i.e. in Eq. (|37|) we set A± 2 = and we 
get the following equation 



In 



T. 



cO 



w 2 



(50) 
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FIG. 5: Contribution of the orbital effect as a function of T c p/T c o for several angles a between H and a;-axis, for r\ — 1.7 (left 
panel), for r\ = 2.5 (middle panel) , and for r\ = 5.1 (right panel). Solid lines are the results obtained for m x /m y = 100 and 
m x /m y = 10; Dashed lines are for m x /m y = 0.01, m x /m y = 0.1. 




FIG. 6: Normalized transition temperature, T c (a) /T c p as a function of a for several T c p/T c o, t/T c o = 0.25, and r\ = 1.7. For 
the purpose of clarity the shown range of T c (a) /T c p is from 0.9 till 1.0. Dashed lines are for A±2 7^ 0. 



Performing average over the Fermi surface results in 

1 



*y/u>* + (Hv F Q) 2 /4 



(51) 



For extremely large magnitude of the external magnetic 
field we can simplify, since u>„ ~ T c q < ujfQ 



In 



t- 



1 



1 



irt" 



1 



T c nT c0 ^ (n + ^) 2 hv F Q 2T c0 HvfQ 



(52) 



Therefore, the upper critical field is (kjv ■ H ^> 



TzhdVF 



From Eq. [52] it is seen that an increase of the external 
field far beyond the value 4>o results in a critical 

temperature T c — > T c q. Hence at high magnetic fields 
the restoration of superconductivity is possible if the de- 
struction of spin-singlet state of Cooper pairs may be 
neglected, as was predicted by Lebedi 2248 Therefore, we 
can infer that within our model the re-entrant phase of 
superconductivity is naturally described. Summarizing 
the above two sections we plot all considered regimes for 
the case of absence of the Zccman effect in Fig. [3] 



H. 



h=0 



2hdv F (T c0 — T c )' 



(53) 
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FIG. 7: The same as in Fig. [5] but calculated for r\ = 2.5, which corresponds to v F = 7.5 x 10 4 m/ sec. 




FIG. 8: The same as in Fig. [5]but calculated for r\ = 5.1, which corresponds to vf = 1.5 x 10 5 m/ sec. 



V. ANISOTROPY OF THE UPPER CRITICAL 
FIELD 

In our numerical investigations we restrict ourselves 
to the following parameters: the inter layer coupling is 
t = 2.27 K, t/T c0 = 0.25, A = 2.8fcT c0 £i and the Fermi 
velocity vp = 5.0 x 10 4 m/ sec<££ Introducing the dimen- 
sionless Fermi velocity parameter, r\ = hvpTrd/4>oHB, this 
value of vp corresponds to rj = 1.7 and d = 1.62 nm42 
The summation over the Matsubara frequencies was per- 
formed numerically. 

Fig. 0] shows the reduced temperature, T c p/T c q, de- 
pendence of the magnetic wave vector HQ^vf /kpTco 



for several values of the Fermi velocity parameter, when 
only the paramagnetic effect is accounted for. Here 
Qc2 = ndH^/fa- The absolute value of the FFLO mod- 
ulation wave vector is also given and it grows from zero 
for T < T* . To highlight the contribution of the orbital 
correction to the superconducting onset temperature, ob- 
tained in the paramagnetic limit, AT c p = T c — T c p, 
and how it depends on the magnitude of the external 
magnetic field applied parallel to the conducting planes 
we performed calculations with Eq. (|45|) . Fig. [5] dis- 
plays the normalized orbital correction, AT c p/T c p, as 
a function of reduced temperature for several angles a 
that the external field makes from the i-axis. The left 
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and middle panels display the results for the velocity 
parameter rj = 1.7 and ?/ = 2.55, respectively. The 
solid lines correspond to the in-plane mass anisotropy 
m x /m y = 100, while the dashed lines display the re- 
sults for m x /m y = 0.01. The right panel illustrates 
the results for rj = 5.1, m x /m y = 10 (solid lines), 
m x /m y = 0.1 (dashed lines). One can distinguish the 
in-plane mass anisotropy from the temperature depen- 
dence of the orbital corrections for angles a ^ ±90°. For 
example, for m x /m y = 100 a decrease of temperature 
from T < 0.9T c o, or an increase of the applied magnetic 
field from H > 0.1-ffpo, first exhibits a weak influence on 
AT cP /T cP , but when T < 0.65T c0 (H > 0.5H P0 . Here 
Hpo = Ao//Ub is the critical magnetic field at T = in 
Pauli limited 2D superconductors) it gradually increases 
|AT c p|, i.e. the orbital suppression of superconductivity 
becomes stronger with magnetic field, when orbital pair- 
breaking is superimposed on the spin pair breaking mech- 
anism. For m x /m y = 0.01 an increase of the applied field 
results first in a progressive increase of |AT c p|. However, 
for T < 0.65X^0 we see an opposite bias, namely strength- 
ening of the applied field rapidly reduces | AT c p\, i.e. the 
orbital pair breaking becomes weaker with the external 
field, and it can almost vanish for some directions of the 
field in the very close vicinity of the tricritical point as 
seen for dashed curves a = 0. For a = 90° the curves 
describing m x /m y = 100 mass anisotropy coincide with 
those for m x /m y = 0.01 and both follow the tendency 
typical for m x /m y = 0.01 mass anisotropy. In Fig. [5] 
both curves are given by the thick lines. We can also 
infer that an increase of the Fermi velocity weakens this 
effect of |AT c p| reduction as seen from the middle panel 
of Fig. El In the FFLO phase, for T < T* , or H > H* , 
the orbital correction in both cases of mass anisotropy 
essentially increases, especially for m x /m y = 0.01 and 
for some angles can show a non-monotonic behavior. The 
further increase of the Fermi velocity can modify the just 
described behavior. Indeed, as seen from the right fig- 
ure the a = 90° curves follow the tendency typical for 
Tn x /m y = 10 mass anisotropy and in the FFLO phase 
they show an upturn. 

Opposite tendency in the field direction dependence of 
the normalized correction, AT c p, for the range of angles 
a = 0°— 70° and for the angles in the close vicinity of 
a = 90° in the case of m x /m y = 100 should result in 
a particular anisotropy of the onset of superconductiv- 
ity. Figs. [6] and [7] show the magnetic field angular de- 
pendence of the normalized superconducting transition 
temperature, T c (a) /T c p 7 calculated at T c p/T c q ~ 0.1, 
0.2, 0.4, 0.54, 0.57, 0.65, 0.84 and 0.99 for the velocity 
parameter rj = 1.7 and rj = 2.55, respectively. In the 
polar plot the direction of each point seen from the ori- 
gin corresponds to the magnetic field direction and the 
distance from the origin corresponds to the normalized 
critical temperature. We see that for m x /m y = 100 the 
reduction of the orbital suppression of superconductiv- 
ity at a = ±90° in the vicinity of the tricritical point is 
accompanied by a grow of cusps at these angles in the 



field-angle dependence of T c (a) /T c p. The cusps appear 
at Q\\Ox, i.e. magnetic field is along the light mass di- 
rection, as intuitively expected, since it is more difficult 
to induce diamagnetic currents with heavier charge carri- 
ers. For m x /m y = 0.01 the overall orbital corrections are 
smaller than that for m x /m y = 100. This is due to the 
fact that in the former case the Fermi surface is smaller 
and hence the diamagnetic response is weaker than that 
in the latter situation. In ^8\T c (a) /T c p is shown for 
ij = 5.1 and m x /m y = 10 (red lines), m x /m y = 0.1 
(green lines). Formation of cusps in the vicinity of the 
tricritical point is also observed, although to a smaller ex- 
tent. In Figs. |6l [7] and [8] the dashed lines are T c (a) /T c p 
obtained for m x /m y = 0.01(0.1 in Fig. [8]) when the r.h.s. 
of Eq. (|37|) is neglected, A±2 = 0. In this case the solu- 
tion (|4"S")) simplifies to 

T c = T cP [1 - At 2 a] (54) 

and such solution is valid for \JtT C Q <C hvpQ, which 
is the beginning of the superconductivity re-entrant 
regime . 22 ' 46 As the charge carrier mass becomes smaller 
the superconducting re-entrant phase begins at a higher 
magnetic field. Since, according to Eq. (|C5j) the sec- 
ond harmonics of the order parameter generates the 
Lowerence-Doniach term in the original expression, Eq. 
(|C9[) . the dashed lines give a hint about its contribution 
to the in-plane anisotropy of the onset of superconduc- 
tivity in layered structures with m x /m y = 0.01 in-plane 
mass anisotropy. We see that the difference between the 
solutions {45]) and ([54]) is negligible for T cP /T cQ ~ 0.57. 
However it is noticeable already for T c p/T c0 ~ 0.65. The 
upper and lower knobs are observed when the full origi- 
nal expression is used, and they are absent for the sim- 
plified version, Eq. (|54[) . So we can infer that the ob- 
served knobs are due to the Lowerence-Doniach term. 
Because this term becomes less important with the field, 
the knobs are absent for T c p/T c q ~ 0.57 and essentially 
pronounced for T c p/T c q ~ 0.85, when m x /m y = 0.01. 
Inversely, for m x /m y = 100 the cusps are profound near 
the tricritical point, insignificant for smaller fields, and 
essentially seen far beyond the tricritical point in the 
FFLO phase. The cusps are induced by the £ 2 a-term, 
which in the conventional phase acquires the following 
form 

t 2 1 
at 2 = ttT cP V — = (55) 

From Fig. ([7]) we can infer that an increase of the Fermi 
velocity leads to a narrowing of the cusp width. However 
such increase of the Fermi velocity makes the cusps less 
pronounced. 

In the FFLO phase hvpQ > T c o, and the solution 
Eq. ([M]) can be used for calculations. The top panels 
of Figs. [6l [7] and [8] illustrate the anisotropy of the super- 
conducting onset temperature in the FFLO phase. We 
see that the cusps induced by the £ 2 a-term becomes even 
more profound with the magnetic field. Moreover, for 
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m x /m y = 0.01 a difference between the results obtained 
within A± 2 7^ and A ±2 = appears. For m x /m y = 0.1 
this discrepancy is also present, although less visible. As 
was shown and explained in Refi^ this deviation this 
time is due to the resonance between FFLO modulation 
wave vector and the interlayer coupling modulated by 
the vector potential. Thus, in addition to the overall 
anisotropy induced by the FFLO modulation and studied 
in Rcfi££, additional cusps develop for certain directions 
of the applied field, when the resonance conditions are 
realized. To describe resonances we have to account for 
the second harmonics, A±2, and then 

s± (Q) s (a + b ± )?+S ± 

In general, in the vicinity of the tricritical point when 
comparing the in-plane anisotropy of T c (a) for the con- 
ventional phase with that in the FFLO modulated phase, 
T < T* or H > H*r& it is obviously seen a significant 
discrepancy. On both sides of the tricritical point, T*, 
the contribution of the i 2 a-term is essential and the ob- 
served difference is purely induced by the appearance of 
the FFLO modulation wave vector. 

The anisotropy of the onset of superconductivity ob- 
tained within our model for t <C hvpQ and m x /m v = 100 
qualitatively similar to that observed in the experiment 
with (TMTSF) 2 C10 4 ^ For H < H* our theoretical cal- 
culations show that in T c (a) /T c p cusps develop along 
the light masses. The same cusps and along the this di- 
rection are visible for H = 20 kOe and H = 25 kOe in the 
experimental data for T c (a) /T c p. Our calculations show 
that for H > H* small dips appear from both sides of 
each cusp. Similar picture is observed in the experiment 
for H > 30 kOe. 

If we compare the field-direction dependence of the 
superconducting onset temperature for T c p/T c q < 0.85, 
valid for HvfQ 3> t, with that in the last panels of Figs. [51 
[7] and [51 where the result of the Ginzburg-Landau regime 
Eq. flTHJ), valid for Hv F Q < t, is shown at T cP /T c0 ~ 0.99 
we see an essential distinction. In the vicinity of T c q the 
anisotropy of the onset of superconductivity shows a typ- 
ical picture for the anisotropic Ginzburg-Landau model. 
T c (a) is maximum for H_LOx near T c o and as seen from 
Fig. [6] also in the vicinity of T* . 



regimes, which we discriminate according to the distinc- 
tive features of the anisotropy of the onset of supercon- 
ductivity and the temperature dependence of the upper 
critical field, (i) In the Ginzburg-Landau regime, when 
H ^ ¥13^7^0, H c2 \ GL ~ (T cP -T c ), the anisotropy 
is well described within the continuous GL model, (ii) 
In the Lowerence-Doniach regime, within tfio/nhdvF <C 
H « VtT^(f> /TThdvF, H c2 \ LD ~ l/V(To-T ct ), the 
anisotropy is mostly determined by the term propor- 
tional to t A /(hvpQ) 2 , which induces knobs in the di- 
rection along the light masses in the field-angle depen- 
deuce of T c ( a), (m) Fo r VtTTo « H « -^-T c0 , 

H C 2\ RS = \J (T c — T c t), the anisotropy is governed by 
the i 2 a-term, which is responsible for the re-entrant of 
superconductivity, (iv) the FFLO phase, H > H* , the 
anisotropy is settled by the interplay between the modu- 
lation and magnetic field wave vectors. The third regime 
can be deep in the four one so the discussed cusps can be 
invisible in the conventional phase. The paramagnetic 
effect is crucial for the description of the upper criti- 
cal field both above and below the tricritical point. If 
the paramagnetic effect is negligible than the extended 
Lowerence-Doniach model restores the re-entrant behav- 
ior with magnetic field originally obtained by Lebed . 22 i 46 
Near T c q the anisotropy of the onset of superconductiv- 
ity shows the smooth variation of T c (a). When reducing 
the temperature, above the tricritical point small cusps 
appear. We may expect that small cusps observed in 
the field-direction dependence of T c (a) /T c p in the ex- 
periment with (TMTSF) 2 C104 near the Pauli limiting 
field, Hpo = 26 kO&i^ could have the re-entrant phase 
origin and are well described by the extended Lawrence- 
Doniach model. A technique that control the anisotropy 
of the upper critical field can provide an invaluable tool 
for investigating the physical origin of the experimentally 
observed upturn of the upper critical field in the low tem- 
perature regime. 
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VI. CONCLUSIONS 



Appendix A: Derivation of the expression for A. 



In this work we have derived the extended Lawrence- 
Doniach model, which allows one to study superconduc- 
tivity of layered materials at high magnetic fields. Within 
this model we have analyzed the field-amplitude and 
the field-direction dependence of the onset of supercon- 
ductivity in layered conductors. Our theoretical analy- 
sis gives rise to the following assertion. There are four 



Substitution of Eq. (JT3J) in Eq. (uDJ results in 



A(r,fc z )ln^- = A(r,fc z ) 

1 cP 

' ttTc ) Y ttT cP 



n MLD A(r,fe), (Al) 
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where we defined a function 



F 



to n (T) 0„ (T) 



(A2) 



When expanding in series, taking into account that 
(T c - T cP ) /T c «lwc obtain 



A ( r , k z) 7F, = A ( r > k z) 



d (h_ 



ttT c p T c 
hn M L D A(r,fc,), (A3) 



T=T C 



and hence 



A (r, k z ) T ° A ^ cP = HmldA (r, k z ) 

where we introduced the following notations 

T c — T c p 



P 



AT r . 



and A is given by 



a-* = i-A *_*.(A 



(A4) 



(A5) 



(A6) 



T=T Cj 



Appendix B: Derivation of Eqs. (|37H39|> 

Solution of the system of coupled equations (|33l - [36]) 
can be found as follows. From Eq. (|36p we find 



/±3 = T- 



tf±2 



L n (±3Q) 
and substituting it into Eq. (|35p gives 

p 



(Bl) 



L n (±2Q) 



/±2±i/±i = A ±2 . (B2) 



L n (±3Q) 

Then substitution of f±\, obtained from Eq. ([34 

, ^ tfp t/ ±2 

J±i = T- — , . „s ±- 



L„(±Q) L„(±Q)' 



(B3) 



when taking into account that within the required ap- 
proximation /o ~ A /L n (q), produces the equation for 
the second harmonic of the pair amplitude, /± 2l 



L n (±2Q) 



P 



L n (±3Q) L„ (±Q) 



L n (0) L„ (±Q) 



/±2 

= A ±2 . (B4) 



Substitution of /±i from Eq. (|B3[) and /±2 ~ 
A± 2 /L n (q± 2Q), obtained within the required approx- 
imation from Eq. (|B4|) . into Eq. (|33f results in the fol- 
lowing equation for fo 



L n (0) + 



P 



P 



L„(+Q) L„(-Q) 

-E 



/o 



Pa 



±2 



L n (±Q)L„ (±2Q) 



A . (B5) 



Since we adopt a second-order approximation in the small 
parameter t/T c0 Eqs. (|B4l - IB5[) acquire the following form 



/o = A 



L n (0) LI (0) L n (+Q) 



L2 (0)i n (-Q) 



E 



PA 



1-2 



f±2 = ^±2 



L n (0)L„ (±Q)L„(±2Q)' 



P 



(B6) 



L„(±2Q) (± 2 Q)L„(±3Q) 



(±2Q)L„(±Q) 



PA r 



L n (0) L„ (±Q) L n (±2Q) ' 



(B7) 



Submitting the obtained expressions for / and /±2 back 
into the self-consistency relation Eq. (Q| results in Eqs. 



Appendix C: Derivation of Eq. ([48]) 

If (hv F Q) < T c0 , or H « ^^^o, then P + P6± « 
<5± and we find from Eq. (|38)) that 



A ±2 « -^A , 



<5± 



(CI) 



with [sec Eqs. (@2J and @J|] 



(C2) 
(C3) 
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where g — hvpQ /Q n . Expansion of these expressions 
with respect to g -C 1 gives 



(hv F QY 



2m 



c± 



2m 



i - 



7(hv F QY 

mi 



and from Eq. (f38|) we find that A±2 reads as 



A 



±2 



{hv F Qy 



■A . 



(C4) 
(C5) 

(C6) 



Substitution of A±2 back into Eq. ([57)1 leads to the fol- 
lowing equation, determining temperature T c of the onset 
of the superconducting state, when the orbital effects of 
the applied magnetic field arc accounted for within the 
second-order approximation in parameter t/T c $, 



P + t 2 a 



(hv F Qy 



(C7) 



where a = 2irT cP J2 n + 9 2 - Making use of the 

expansion of a into a series 



fi?, 



i {hv F QY 

8 SI* 



we obtain equation for T c 



1 _i(hv F QY 



(C8) 



sil {hv F QY 



(C9) 

After introducing T ct , as it is done in Ref. (f2"r?)) , which 
accounts for the coupling between adjacent layers, finally 
we obtain Eq. 
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